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I . SUMMARY 


A set of general equations for evaluating the sensitivity of the 
differential Absorption and d^attering (DAS) technique based upon a conventional 
analysis of statistical errors is derived in this report. The equations are 
put in a proper form for evaluating total column density and range resolved 
concentration measurements of a variety of atmospheric species. The derived 
equations are subsequently used to analyze the sensitivity of DAS in three 
specific applications assuming realistic parameters for the optical and elec- 
tronic components of proposed DAS systems. 

The three DAS applications evaluated are: (1) measurement of nitrogen 

dioxide (NO2) at ground levels over a horizontal path, (2) measurement of atmos- 
pheric ozone (O3) depletion in the wake of a jet engine at 20 km altitude, and 
(3) measurements of the ozone distribution in the atmosphere from an orbiting 
space platform, in a downward viewing mode. The results of this study have 
shown that with reasonable laser energy and telescope receiver dimensions, DAS 
is capable of meeting requirements for performing these measurements. 

II. INTRODUCTION 

IDifferential Absorption and Scattering (DAS) is a single ended active 
remote sensing technique for monitoring concentration of various atmospheric 
gases. Since DAS relies on Rayleigh and Mie backscatter, ranging information 
is possible. The technique was first applied by Schotland (1965)^ to determine 
the vertical profile of atmospheric water vapor using a ground based ruby 
laser which was tuned through a water vapor absorption line. Although a detailed 
uncertainty analysis was not performed for this early measurement, a signal-to- 


noise ratio for the return signals was estimated and standard deviations were 
obtained for sets of measurements , The results of this early measurement did 
show a trend of increasing cumulative water vapor with altitude as expected. 

A recent report of Measures (1971)^ illustrates the advantages of differential 
absorption to either Raman Scattering or Resonant Fluorescence. This result 
is expanded in a later report by Measures and Pilon (1972]^ but no detailed 
sensitivity expression was derived. 

Kildal and Byer (1971)^ presented details of comparison for Raman back- 
scatter, resonance backscatter, and resonance absorption. The resonance 
absorption considered requires a remote detector on a retroref lector and does 
not include the phenomena of molecular or aerosol backscatter utilized in DAS 
measurements . The DAS concept was included in a later report of Byer and 
Garbuny (1973),^ Other aspects of an error analysis of the DAS, (or DASE for 
differential Absorption of Scattered d^ergy) are presented by Schotland (1974)® 
and Ahmed (1973).^ In the work presented here a more basic approach was used 
than in the earlier works and the analysis is developed around the concept of 
a system that measures energy (or photons) rather than power with a given 
bandwidth. The measurability criteria or sensitivity used in this report is a 
predicted standard deviation based upon Gauss and Poisson statistics. A per- 
centage predicted uncertainty is ejq)ressed by dividing the propagated standard 
deviation by the ’’measured" value. There is a 67% probability that a single 
measurement of a given quantity (defined as the* signal S in signal-to-noise S/N) 
will fall within one standard deviation (the noise N) of its true value. The 
expressions derived in this report have been incorporated into a DAS simulation 
computer program developed by Stanford Research Institute in order to more 
justifiably predict the retrievability of pollutant profiles in the atmosphere. 
The equations derived allow one to determine the sensitivity of measuring a 
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gas concentration as determined by limitations of background noise, attenuation 
due to factors other than the gas of interest and backs cattering variations. 


III. THE DIFFERENTIAL ABSORPTION AND SCATTERING TECHNIQUE 


The concept of differential absorption and scattering is illustrated in 

the case treated by Byer and Garbuny^ shown here as Figure 1 through Figure 3. 

A laser beam is encompassed by the field of view of the telescope as shown in 

Figure 1. The cloud at the center of Figure 1 represents a gaseous plume for 

which a concentration measurement is desired. Two range cells corresponding to 

detector gate time intervals T. and T. , at ranges R. and R contain 

1 1+1 1 3+1 

Rayleigh and Mie scatterers such that backscattered laser light plus background 
(e.g., scattered sunlight) are received at the detector situated at the focal 
point of the telescope. The AR, shown in Figure 1 represents the spatial 
range over which a difference in column content applies in determining concen- 
tration as a function of range. The expected return signal is plotted in Figure 
2 as a function of range for the strongly absorbed signal Sj and the weakly 
absorbed signal S^. Note that the R^ dependence has been removed from the 
return signal profiles. 

The logarithm of the return signals Logg[Si/S2] at each wavelength is 
plotted as M in Figure 3 illustrating the column content as a function of 
range under conditions of identical losses at both wavelengths due to materials 
other than the particular gas of interest. 

Finally, the curve labeled as N in Figure 3 depicts the constituent (in 
this case a pollutant) gas concentration as a function of range as obtained by 
differentiating M with respect to range R. The precision with which concen- 
tration can be determined is a function of the noise on the signal return and 


the range resolution or averaging time 


T. 

3 


required. 
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A. The Lidar Equation 

In this section an expression for the backscattered signal is defined for 
a DAS system in terms of atmospheric and system parameters. The Lidar return 
signal strength S. . observed at the detector from range cell j at range R 

3 

for wavelength X. is given in photons per range cell (for Q = 1) by* 

1 

QUi A cT^ ( ) 

" h^ ^ 2 ~ ^5 N(r)a. (r)] dr| + Qn^p + Qn^g 

Cl) 

where 

Q represents electronic amplification which converts photons received 

per range cell (T^ ) to charge stored by electronics; 
is the energy of the laser pulse at wavelength joules; 

h is Planck's constant, joule-sec; 

A is the effective area of the receiver telescope, m^; 

Rj is the range at which the j^^ range cell is located, m; 

is the combined Rayleigh and Mie backscattering coefficient at 
wavelength X^^ in range cell j, m“^ sr~^; 


* It is assumed that 



where 


a 


"j 



and ' is the backscattered laser signal at wavelength X^^ from range r 

observed at the lidar receiver. This assumption requires that the average return 
signal strength from range cell j can be represented by the return signal 
strength calculated at the midpoint of the cell. 
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c is the speed of light, m/sec; 

Tj is the integration or gate time for the detector, corresponding 

to a range cell at Rj , sec; 

Tg is the optical efficiency of the telescope-filter-detector system 

at wavelength 

C^(r) is the extinction due to interferents such as Rayleigh scattering, 
Mie scattering and atmospheric constituents other than those 
included in N(r)a^Cr) at wavelength and range r, m*"^; 

NCr) is the concentration of atmospheric constituent gas to be measured, 
atm; 

Uj^Cr) is the photo-absorption coefficient at wavelength X^, (atm cm)^^ 

n^j^ is the equivalent number of photons counted due to the detector dark 

current during time , photons; 

n^g is the equivalent number of photons counted due to background such 
as scattered sunlight, photons. 

The return signal strength due only to backscattered laser light is obtained 
by subtracting the background and detector noise levels defined as 

Si = QCnip + Hig) (2) 

from the return signal strength of Equation (1) . The Lidar equation evolves by 
subtracting Equation (2) from Equation [1) and dividing by the electronic ampli- 
fication factor Q to obtain 

n^j = EXP |-2J^^[?i(r) +NCr)a^(r)] dr | (3) 

where n is photons per range cell resulting only from backscattered laser 

ij 

photons . 
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B. The Column Content Measurement 


The column content of a given gas is defined here as the integral of gas 
concentration over a specified distance. For example, a molecular pollutant 
having a uniformly distributed density of N atmospheres 

(1 atmosphere = 2.69 x 10^^ molecules / cm^) has a column content over a range 

of 1 meter given by 100 N atm cm or over a range of 1 km given by 10^ N atm cm. 

The column content. of a gas to range can be obtained. by measuring 

n, . at Ai (on the absorption peak) and n^ . at X2 (off the absorption peak) 

Ij 

and solving the ratio of n^^^ to n^^ of Equation (3) for the integrated 
concentration (assuming a constant difference in absorption coefficient) to 
obtain : 


M. =/^N 


(r)dr 


2 (a 


1 ( n- .U, 3 1 -T ,An J?i ) 

■ ' 1 {Log - 2 lA * Logg - 2/ [ 5 iCr) - 52 W] dr[ . 

1 - ^ P2j^s2'^2 o ) 

C4) 


The dominant term of this calculation must be Log^(n2j U2) for any 

measurement so as to minimize the requirement of a precise a priori knowledge 
of the other terms. Requirements on precision and resulting limitations are 
discussed in a following section. The column content as expressed by Equation (4) 
is valid also when a retroref lector is used where replaced by the 

ratio of reflectivity at the two wavelengths. 

In terms of the measured parameters , S^* and Equation (4) becomes; 

(So-; - S2) Ui ^l^li^sl ^ 

Loge yj + Log - 2/ ci(r) - S2(r) dr}. 

CS^j - Si) U2 ® ^2823^52 *'0 ) 


“ 2(ai 


- «2) .{' 


C5) 


In developing equations (4) and (5) it is assumed that photo -absorption has no 
r dependence which is not necessarily true such as when atmospheric pressure 
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or temperature changes with r. When this approximation is not valid the proper- 
ties of the measurement path must dictate how oj and 02 will be determined as 

weighted averages. An alternative would be to express "pollutant" level as an 

R 

optical depth difference (i.e., f N(r)[ai(r) - a 2 (r)]dr) which is readily 

o 

obtained from Equation (3) . 


C. The Range Resolved Concentration Measurement 

The concentration of an atmosphere constituent can be determined as a 
function of range by dividing the difference in column content at two ranges by 
the difference in the respective ranges. The resulting measure of concentration 
will represent a mean over the range increment involved and will have a range 
position determined as the mean range. 

The average concentration <N>^ contained in a range increment R^ to ^j+i 
is defined as 


<N> 


I^+l 


Then from Equation (5) it follows that 


( 6 ) 


<N> 


where 


1 


2AaAR- 

J 


(Sij^.1 - SOCSj. - S2)J^ 



Logs 


Ij 




Aa = Oj - «£ 

and 


(7) 


AR. = R. , - R 
3 3+1 j 


The values of term © are of primary importance 
since they are the largest determining factor of 


for an uncertainty analysis 
Equation (7) whereas the values 
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in will exactly cancel when Rayleigh scattering dominates or when Mie 
scattering dominates in both range cells j and j+1. The values in the © 
term represent other loss mechanisms and must be Included for completeness. 
Limitations due to terms © and © are discussed below. 


IV. PREDICTIONS OF DAS SENSITIVITY 


An estimate of uncertainty, therefore sensitivity, of a DAS measurement 
can be obtained by use of the principle of least squares. Error propagation by 
least squares is treated by Bevington [1969)® and other standard texts on 
statistics and error analysis. The least square procedure is described briefly 
below. Assume a normal distribution of error 5 (the standard deviation) 
about the true value of a set of n measured parameters describing a function 

f (X 2 .... • Consider a Taylor expansion of f about the mean values of 

each X^; square f - f^j ayerage over a large number N of observations to 
eliminate odd order terms, then drop fourth order and higher terms to obtain 
the error propagation form of least squares 

[«fCXi . . . [(!-)' (6X^2] C8) 

where 

6f(Xi . . . X^) = 

and 

5X. = 

1 

The zero subscript denotes mean value. The above definitions of 5f and 5 X^ 
show their relation to variance and standard deviation. For purposes of 
uncertainty prediction the values of 6 Xj^ can be determined from Poisson 




8 



statistics, from standard deviations of measured data, or from other statistical 
arguments that fit a given situation. 


A. Predicted Sensitivity in Column Content 


The predicted uncertainty in column content is obtained by operating 

on Equation (4) with Equation ^8} to obtain 

2 


5M. = U- 

J I 2Aa CS 2 j " 


C<SS2)^ _ (5S 

^Si 


lj)2 + (SSi)^ ^ 

3ij - Si)2 \U7/ 


(5 


4(6 (CR) 




C<S Ag) ^ 2 

(Aa)2 


1/2 


(9) 


The values of Tgj and Tg 2 are not subject to random fluctuations and are 
not included in Equation (9). The values of 6S^, 6(3i^j/32j) and 6 (CR)j 

are discussed below. The value of 5 (Aa) is considered inconsequential for 
applications considered in this report but for other applications standard devia- 
tions of measured values may be required. The uncertainties in laser energy 6 
can be made inconsequential by design but are included in Equation (9j to indicate 
the precision with which must be measured. 


B. Predicted Sensitivity in Range Resolved Concentration 

The predicted uncertainty in range resolved concentration measurement 5 Nj 
is obtained by operating on Equation (7) with Equation (8) where f = <N>j to 
obtain 



The value of 6 (Aa) is again ignored in this report for the same reasons stated 
in the above section. Suitable assignments for the rest of the standa:?d devia^ 
tion terms will now be discussed. 
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C. Uncertainty in return signal measurements 

The uncertainties in S^j and are obtained by operating with Equation 

(8) on and as defined by Equations (1), (2), and (3) to obtain 

+ C5nig)2] + (SQ)^ CniR * "ID + "Ib)^ dD 


and 


|SSi)2 = Q2[(SniD)^ * PniB)^] + C«Q)^ (nio + n^B)^ 


(12) 


The 5 ’s of Equations (9) » (10) » (11)> and (12) represent standard devia- 
tions from the mean. The standard deviations in laser backscatter Snj^j. is 
assumed to be Poisson such that <S n^j^ = /n^j^ and can be calculated for an 
assumed laser lidar system and gas constituent profile from Equation (4) . The 
value of the noise resulting from the detector varies from 0 (with 

respect to one photon per range cell) for a cooled photomultiplier to a value 
determined by^ 




/Tj 

D*^h^V-: 


(13) 


for a solid state detector. The noise equivalent power (NEP) is a parameter 
determined by the selected detector. The value of NEP depends upon the active 
surface area A, time response Tj and detectivity D* of the selected 
detector. Typical values of D* are on the orderly of 10^^ w~l cm sec”"^/^ 
whereas the theoretical limit can be as high as 10^^ at 3.0 u (see, for 
instance, Putley^). 

« 

The noise generated from background fluctuations are assumed to be 
related to the background n^g by Poisson statistics. For a system opto- 
electronic efficiency n, telescope field of view 0^, telescope receiver area 
A, detector gate time Tj , spectral bandpass AX^, and a spectral radiance 
in w m“^ sr“^ nm”l the received photon flux from the sky background is^ 


(5niB)2 


^iB 


iTJT9^ATjAXiBr 

hvi 


(14) 
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Values o£ Bj. for the region 300.0 nm to 420.0 nm have been reported by 
Knestrick and Curcio^^ and range up to 9 10“^ w sr“^ nm“^ . 

Uncertainties introduced by electronic amplification factor Q will depend 
upon the selected system. Specific knowledge of Q is not required since all 
terms containing it can be expressed in terms of 6Q/Q. Typical values of 6 Q/Q 
can range from <1% for a gated integrator coupled to an analog-to-digital con- 
verter to 2.3% for a logarithmic amplifier. 

D. Uncertainty in Backscatter 

Backscatter results from a sum of the Rayleigh and Mie backscatter which 

vary with and (1 ^ ot ^ 4) , respectively, ignoring resonances. The 

value of Bjj/3 2 j becomes a limiting factor in the calculation of column content 

Mj if the absolute value of 3 Logg(Xi/A2) is on the order of 2MjAa, Laser 

wavelengths should be chosen carefully so as to maximize 2M.Aa - 3jLog CX 1 /X 2 ) | 

J e 

while maintaining 6M^ of Equation (9) within acceptable limits. It is possible 
to correct for this difference in backscattering between wavelengths to some 
degree of precision as indicated by the form of Equation C4) . The value of 

niust now reflect the precision to which the backscattering is known 
to be non-Rayleigh when a Rayleigh scattering relationship is assumed. The amount 
of Mie scattering can be estimated from visibility by the technique of Fenn.^^ 

When a retroreflector is used the value of 5 C5ij/62j) must reflect the 
precision to which the wavelength dependence of the retroreflector is known. 

For range resolved measurements the limitations due to backscattering are 
in general less severe due to the self-cancelling nature as seen by inspection 
of Equation (7) . The value of the © term is exactly zero if the backscatter 
wavelength dependence is unchanged between range cells j and j+1. This 
condition will exist if scattering is dominated by Rayleigh scattering or by 
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Mie scattering Cwith identical wavelength dependence) in both range cells. 

When backscattering o£ Rayleigh and Mie scattering or the wavelength dependence 
of Mie scattering (for the Mie dominant case) varies between range cells j 
and j+1 the value of term © of Equation (7) will be non-zero. Under this 
latter condition the value of 53j'/3j' must reflect the precision to which 
C6i j +1 j+1 ^Ij^ known. Under the conditions for which term 0 

of Equation (7) is zero the value of 5Bj'/Sj’ must reflect the precision to 
which identical wavelength dependence of scattering in range cell j and 
range cell j+1 is known. 

E. Uncertainty in Extinction 

The column content measurement as obtained via Equation (5) is limited 
when extinction due to secondary mechanisms such as absorption and scattering 
by other gases becomes comparable to the term containing the return signals 
The effects of these other mechanisms are reflected in the Ci(r) - C 2 (r) term 
in Equation (5) and can be systematically corrected when known. When this 
systematic correction (ClCr) - CzC^)) is required the measurement limitation 
becomes dependent upon the precision with which this correction can be determined 
as reflected in the 6 CCR)j term of Equation (9), In the case of column con- 
tent measurements the correction applies to the entire column, whereas for the 
range resolved concentration measurement it applies only in the range increment 
ARj in which the concentration is being determined. 

The required systematic correction for secondary loss mechanisms in range 
resolved concentration measurements is contained in Equation (7) term 0 and 
becomes important when term is comparable to term containing the return 
signal strengths. The precision to which Si(r) - CaCt) is determinable then 
becomes the limiting factor in measuring <N>^ as reflected by 6<?R>^ in 
Equation (10) . 
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If secondary losses result primarily from another gas the value of 6<CR>j is 
expanded in terms of the gas concentration and cross section uncertainties. If 
the limitation is due to aerosol or Rayleigh scattering the 5<?R>j must be 
evaluated in terms of the precision to which the scattering is known. For 
purposes of modeling, estimates of atmosphere constituency can be obtained from 
Elterman, ^ McCormick, McClatchey, et al,,^^ and NASA publication SP-285.^^ 
The results presented in the following section evaluate feasibility in 
terms of backs cattered signal strength only in a pure Rayleigh atmosphere and 
do not consider backscatter limitations or extinction from secondary mechanisms. 
These considerations should be included as the next degree of sophistication. 

V. RESULTS 

A. Nitrogen Dioxide Measurements at Ground Level 

The configuration used to determine predicted standard deviations for DAS 
measurements of NO 2 is the same as that described in Section III. In this 
section a parametric study is included, the object of which is to develop an 
intuitive concept for the impact of various experimental parameters on the 
measurement sensitivity. The absorption coefficients used are 17.9 (atm cm) 

at 448.0 nm and 10.5 (atm cm) at 446.0 nm as obtained from the resolved 

0 

absorption spectra of Hall and Blacet^® shown in Figure 4. A higher resolution 
absorption coefficient measurement could markedly increase the ranges of NO 2 
concentrations considered as measurable by increasing Aa of Equation (9) . The 
total Rayleigh-Mie extinction coefficient is taken as 5 = 1.0 ^ 10”*^ m“^ 
whereas the backscatter coefficient is assumed as 3 = 4,0 x 10“^ m“^ sr”^ as 
suggested by Measures,^ The detection system is assumed to have a 6-inch 
diameter telescope, an overall quantum efficiency of = .20 and a detector 
gate window (or reciprocal bandwidth) of T. = 10 ns. The laser power is assumed 
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to be 75 kW with a pulse width o£ 10 ns yielding an energy per pulse of 


ui = U 2 = .75 millijoule. A background intensity of n^^g = n 2 g = 100 photons^ 
per gate time is assumed for these calculations to represent daytime measurements. 


The above set of parameters is not intended to represent an optimized system but 
to provide a base from which an optimized system can be defined as described 


later. 


The predicted standard deviation of an average concentration of NO 2 over 
range R is shown in Figure 5 for various ranges as a function of NO 2 concen- 
tration for the above conditions. From Figure 5 it is apparent that from a 
10 ns range cell at 900^m the column content can be measured to no better than 
about 75% and then only at an average concentration of 1.1 ppm over the 
900 m. The degree of sensitivity Cpredicted percent standard deviation) 
improvement is shown on Figure 5 as the distance from the LIDAR to the range 
cell decreases to 100 m showing a broadening of the measurable set of concentra- 
tions. At 100 m the return signals from the range cell are capable of determining 
the average NO 2 concentration to within 20% if it falls in the region of 1 to 
45 ppm. In practice upper or lower limits can be set on an off-scale measure- 
ment by determining whether the predicted error is large due to strongly attenuated 
signals or a lack of signal differential, respectively. The predicted standard 
deviations are expected to improve for a night operation where the background 
would be zero. We can define a span of range and concentration within which 
meaningful measurements of total burden can be perfcyrmed by setting Equation (9) 
equal to the required minimum of precision [say, AM/M = 20%) in column content 
and solving for M at selected values of range R. 

The regions of acceptable NO 2 concentration measurements (AM/M < 20%) are 
shown in Figure 6. The region enclosed by the curve A bounded by 1.1 ppm, 

45 ppm, 100 m, and 400 m is that region in which NO 2 concentration can be meas- 
ured to within 20% accuracy or less using the set of conditions assumed earlier. 
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The increase in the region of acceptable accuracy caused by increasing the 
absorption at the peak and decreasing the absorption at the minimum by a factor 
of five is shown by curve D bounded by 11 ppm, 0.7 ppm, 100 m, and 1000 m. 
Another means of enhancing the range of detectable concentrations is to '^extend 
the integration (or detection gate) time. The results of enhancement of this 
nature by an order of magnitude is illustrated by the increased region of 
detectability outlined by curve B, Further improvements can be obtained by 
increasing the laser power by an order of magnitude as shown by curve C. 
Improvement of curve C over Curve B is caused by the background which 
remains constant with a power increase but increases with increased gate time 
or number of pulses , 

Now we can consider viable alternatives to improve the design of the 
system. Increasing power as opposed to increasing gate time (or averaging over 
a large number of pulses) is advantageous because of background interaction. 

There is also support for seeking a narrow linewidth laser if the more resolved 
structure produces increased absorption coefficient differences. In addition, 
narrow linewidth filters compatible with a narrower laser linewidth will decrease 
background when background is a limiting factor. However, the laser output energy 
per shot decreases as linewidth is narrowed indicating that an optimum combination 
is reached when laser energy losses due to narrowing become greater than the 
increase in absorption coefficient difference (Act) . 

B. Ozone Concentration in a Jet Wake at 20 km Altitude 

In this section the DAS technique is analyzed for feasibility of detecting 
a depleted region of O 3 due to a catalytic reaction by NO^ in the wake of 
a jet engine. The catalytic reduction of O 3 by NO^, identified by Johnston^^ 
in 1971, can be a detrimental result of SST operation in the stratosphere. 
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The consequence of such an occurrence would be an increase in ultraviolet sun- 
light intensity at the earth’s surface since ozone is a very effective ultraviolet 
filter. 

Consider a remote laser DAS system operating from an aircraft or balloon 
as illustrated in Figure 7, Assume an ambient O 3 concentration of 
16 X 10“® atm or 2 ppm (_at 20 km the atmospheric density is .079 atm), a 
Rayleigh scattering coefficient of .01 km“*^ at 300.0 nm (Figure 8 ) and 
5.57 X 10 ^ km“^ at 347.2 nm (Figure 9), a backscatter factor^® of .12 sr“^ , 
a detection efficiency of ,05, detector gate time of .16 ysec (25 m range cells) 
and a telescope diameter of .25 m. The scattering coefficients given in Figures 
8 and 9 are discussed in more detail in the following subsection. The lasers 
are assumed to have 1 joule per pulse at 300.0 nm and 347.2 nm for which O 3 
absorption coefficients are 10. 8 (atm cm)^“^ and .014 (atm cm)^“^, respectively. 
The jet wake is assumed to be 100 m wide located at a range of 500 m from the 
observation platform as shown in Figure 7. Substituting the above parameters 
into Equation (3) the number of photons received at the detector from each range 
cell are calculated and listed in columns 3 and 4 of Table 1. The uncertainty 
in number density of O 3 for 25 m and 100 m range cells are listed in columns 
5 and 6 respectively. The' results show that O 3 concentrations of 5.4 x 10“® atm 
to 8.9 X 10~® atm can be detected using the 25 m range cells. The precision 
can be increased to 8.0 x 10 ”^ to 1.0 x 10 ~® as noted in column 6 by using 
100 m range cells. Although we increase precision by using 100 m range cells 
there is a corresponding loss in spatial resolution of locating concentration 
changes. Inclusion of particulates in the model would not significantly improve 
feasibility since it is not weak signal limited. 

The conclusions to be derived from Table 1 then are that in monitoring O 3 
concentration in a jet wake a decrease in the O 3 concentration by factors 3 to 
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10 can be observed under the simulated conditions. It should also be noted 


however, that increasing the absorption cross section by using shorter wave- 
lengths [(e.g., 0 “ 300 (atm cm) at 250.0 nm] would also increase the 

precision of such a measurement. A considerable decrease in return signal 
n^^^ can be tolerated without increasing the noise appreciably since is 

so large as shown in Table 1. In a recent calculation Brockman and Seals^^ 
determined that a two ended absorption measurement in the infrared with a tunable 
diode laser provides sensitivities of .02 to .2 ppm over 10 m for O3, CO, NO, 

H2O, and SO2 with minimal interference effects. This result corresponds to a 
sensitivity of 1 to 10% of ambient for 10 m resolution which is more sensitive 
than the sensitivity of 50% to 6% calculated here for 25 m to 100 m resolution, 
respectively. The DAS technique does allow the measurement to be performed 
remotely rather than flying through the wake as required by a two-ended absorp- 
tion measurement and DAS provides ranging information. 


C. Ozone Measurements from an Orbiting Shuttle 

The schematic arrangement considered for ozone profile measurements from 
an orbiting shuttle at 185 km altitude is shown in Figure 10. The 15-km deep 
ozone peak above 15 km in Figure 10 corresponds to the peak of the ozone concen- 
tration shown in Figure 11 as determined for the Midlatitude Summer Model 
Atmosphere. 

A scattering length of 4 km corresponding to a gate time of 26.67 micro- 
seconds is applied to the observed photomultiplier signal. The 4-km deep 
scattering lengths vary from 5 km to 45 km altitude at 4-km intervals. The ozone 
absorption cross sections used for the calculations shown in Figure 13 were 
measured by Griggs^^ with a resolution of .2 nm. 
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To perforin sample calculations only Rayleigh scattering as obtained from 
"Optical Properties of the Atmosphere (Revised) is used for the backscattering 
medium of this DAS simulation. Consider the wavelengths of 300.0 nm (doubled 
dye laser and 347.2 nm (doubled ruby laser) as the laser wavelengths for 
strong 10.8 (atm cm)^'^ and weak 0.014 (atm cm)^"^ O3 absorption, respectively. 
Using the Rayleigh scattering data for Midlatitude Summer at 337.1 nm and 
correcting with the scattering coefficients and integrated scattering 

coefficient at 300.0 nm and 347.2 nm are calculated and plotted in Figure 8 and 
Figure 9, respectively. The backscatter as obtained from McClatchey, et al.^^ 

is .12 sr~^ of the scattering coefficient at range R (altitude h) . Each 

laser pulse is considered to have 1 joule of energy. The combined efficiency of 
the optics and photomultiplier is considered to be 5%. The diameter of the 
receiver telescope is assumed as 1 m. The uncertainties reported here do not 
include multiple pulse averaging. Using the conditions defined above and values 

of / (r) dr from Figure 9 and Figure 10 the predicted return signal is 

*^o 

calculated with Equation (4) at 4-km intervals for each wavelength as listed in 
Table 2. 

The predicted standard deviations (uncertainties) on column content are 
calculated using Equation (9) and listed in Table 2. These results are also 1 
plotted on Figure 13 as error bars on the assumed distribution of column content 
as seen from the shuttle to various range cell al'titudes. Note that the 
predicted uncertainty, which is a measure of predicted sensitivity comparable 
to the inverse of signal-to-noise, increases from 3.0 m atm cm at 42 km to 
8.0 m atm cm at 6 km. The percentage standard deviation is greater than 100% at 
46 km and decreases to a minimum value of 2% at 22 km then increases to 3% at 6 km 
altitude. The percentage standard deviation is dominated by the distribution 
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of the column content M being measured; however, for evaluating techniques of 
measuring ozone content of the atmosphere it is the more important term to 
evaluate. 

The standard deviations in ozone concentration are evaluated with Equation 
(10) and listed in Table 2 as AN and plotted as predicted error bars on 
Figure 11. The result yields an uncertainty of 13 (relative to 13) atm 
at 44 km (or 100%) then diminishes to 7 (out of 89) atm (or 8%) 
at 32 km then increases again to 25 (out of 53) atm (or 47%) at 
12 km altitude. From the signal return strengths n^^j and n 2 j of Table 2 
it is apparent that the precision limitation in the upper atmosphere is limited 
by the amount of backscattering obtained whereas at lower altitudes attenuation 
at the strongly absorbed wavelength drives up the standard deviation. It is 
apparent from the n 2 j return signal that the loss is more than compensated 

for by the increased backscatter caused by increasing air density. 

The solutions presented assume no background fluctuations (Sn.^ = 0), no 

iB 

detector noise (5 n^^ = 0) and no amplifier noise (5 Q/Q = 0), These preliminary 
calculations were performed on a programmable desk top calculator with limited 
capacity. The equations are being coded onto a large computer subsequent to 
which a variety of such problems will be analyzed for feasibility including all 
uncertainty estimates. Indications are that the noise terms included here are 
the dominant ones. ' i 

Results of the analysis at this point support the feasibility of monitoring 
the ozone profile from a shuttle. Some situations under consideration for 
future analysis include double peaked ozone distributions, a terminated signal 
return due to clouds, and total column contents of such gases as SO 2 and NO 2 . 


VI. CONCLUSIONS 


The analysis technique presented in this report is capable of predicting 
the sensitivity attainable with a specified lidar system for a given atmospheric 
model using Differential Absorption and Scattering, Such calculations are 
necessary for determining feasibility of developing the technique for a given 
application and for estimating design parameters of a DAS Lidar system. Capa- 
bilities of the DAS technique for three applications have been analyzed and 
can serve as initial design criteria. 

The analysis of horizontal NO 2 concentration measurements indicates that 
the proposed system can measure average NO 2 concentrations (or column contents) 
to 400 m over a concentration range of 1.1 to 40 ppm to within 20%, which is in 
the region of interest of the recommended threshold level of 5 ppm,^^ The 
design parameters analyzed in this report indicate that in order to monitor 
ambient NO 2 concentrations in a range of .01 to .2 ppb^^ lasers of considerably 
more energy than .75 mj per pulse and/or integration times of much greater than 
10 ns are required. In addition, the return signals can be interpreted to 
indicate whether an out of range measurement (>20% predicted standard deviation) 
is above or below the concentrations allowed for reasonable errors. 

The analysis of ozone concentration in a jet wake at 20 km indicates that 
a drop to about 30% of ambient can be detected which is sufficient to indicate 
the presence of the wake and thus determine its persistence and spread with 
time. The analysis of an ozone profile measurement from shuttle at 185 km has 
shown that for the system described we can expect to determine ozone concentra- 
tion to within 10% at the peak (h = 22 km) and 68% at 8 km and 41 km with con- 
centrations at altitudes between 14 and 40 km being measurable to within a 
precision between 6.3% and 36%, The same system analyzed for total burden to 
the satellite measurements is predicted to have a precision range of 2% at 17 km 
to 11% at 38 km with all precisions below 38 km being less than 11%, 
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The next step is to include interference effects from other molecular gases 
that have absorption structure within the wavelengths used for the measurement 
based on an approximate model of gas distribution over the intended absorption 
path of the measurement. The analysis has been based on a photon counting 
system and neglects statistical errors caused by amplifier gain when large 
count rates are encountered. The calculations are being expanded to include 
these sources of error. In a similar fashion, this evaluation technique can 
be extended to infrared systems where photon counting is not possible. 
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Table 1. Ozone Density, Estimated Signal Returns, and Predicted 
Uncertainty for Ozone Concentrations as a Function of 
Range for Jet Wake Concentration Measurements. 


Range 

(m) 

O 3 Density 
( 10 "® atm) 

Ni (3000) 
(counts) 

ANq 

N 2 (3472) qq- 9 atjn) 

(counts) AR = 25 m 

ANq 

( 10 "® atm) 
AR = 100 m 

400 

16 

576076 

437790 \ 







J 5.4 




425 

16 

505653 

387686 ) 




450 

16 

446926 

345704 




475 

16 

397470 

310179 








N 


) .81 

500 

8 

356992 

279856 




525 

0 

323640 

253765 




550 

0 

294740 

231154 ) 







> 7.1 




575 

0 

269532 

211430 j 


/ 







. .98 

600 

8 

246349 

194122 




625 

16 

224970 

178850 




650 

16 

206105 

165308 ) 







> 8.9 




675 

16 

189382 

153244 ) 

> 




Other Pertinent Parameters 


n = .05 
E = 1 

T = .16 X 10"® sec 
d = .25 m 

Xj = 3000 X 10**10 m 

01 = 1.0 X 10"® m"l X .12 


Cj = 1.0 X 10“5 m"^ 

01 = 10.8 atm"^ cm~^ 

X 2 = 3472 X 10~10 m 

02 = 5.7 X 10“® m'^ X .12 
C2 = 5.7 X 10"® m“l 

02 = .014 atm"^ cm“^ 
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Table 2. Calculated values of received signal strength on resonance absorption (nij), received 
signal strength off resonance absorption (n 2 j), predicted standard deviation in ozone 
column content (5M) and predicted standard deviation in ozone concentration (6N) as 
functions of range (R) from the shuttle. 


Altitude 

h(km) 

Range 
Rj (km) 

nij 

(Photons per 
range cell) 

^2j 

(Photons per 
range cell) 

M 

(matm cm) 

5M 

(matm cm) 

<N>j 
(n atm) 

6N 

(n atm) 

6 

179 

22 

18829 

316 

9 



8 

177 





47 

31 

10 

175 

36 

17852 

297 

7 



12 

173 





53 

25 

14 

171 

48 

13515 

276 

6 



16 

169 





76 

22 

18 

167 

62 

8532 

245 

5 



20 

165 





128 

18 

22 

163 

115 

5095 

194 

4 



24 

161 





169 

12 

26 

159 

299 

3029 

126 

2 



28 

157 





128 

8 

30 

155 

555 

1847 

75 

2 



32 

153 





89 

7 

34 

151 

695 

1064 

39 

2 



36 

149 





54 

8 

38 

147 

641 

613 

18 

2 



40 

145 





28 

10 

42 

143 

470 

352 

6 

3 



44 

141 





13 

13 

46 

139 

302 

202 

1 

4 






Figure 1. Simplified diagram of the Differential Absorption and Scattering (DAS) LIDAR Technique. 
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Figure 2. Simulated two-wavelength signal strength received from a DAS LIDAR system as described 
in Figure 1. 
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Figure 3. Simulated results from processed DAS LIDAR return signals shovm in Figure 2. 
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PARAMETRIC STUDY OF 
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Figure 6. Parametric study of system constants to evaluate techniques of en- 
hancement for NO 2 concentration measurements by DAS. 
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AIRCRAFT JET WAKE 

MONITOR 
SYSTEM 


• NO CATALYTIC DEPLETION OF 0 

• 22 km ALTITUDE 

• 25 m RESOLUTION 


Figure 7. System schematic to monitor ozone depletion at SST altitudes from 
an aircraft using DAS LIDAR. 
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Figure 8. Rayleigh extinction coefficients and integrated extinction coef- 
ficients at 300.0 nm plotted as a function of altitude but observed 
from the NASA shuttle. 
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RAYLEIGH EXTINCTION COEFF. (km ) OR INTEGRATED EXTINCTION COEFF, 


1 


RAYLEIGH EXTINCTION 
AT 347.2 nmTHROUGH THE 
ATMOSPHERE FROM 185 km. 


.185 


■ 0 10 20 30 40 50 

ALTITUDE h (km) 

Figure 9. Rayleigh extinction coefficients and integrated extinction coef- 
ficients at 3472 ran plotted as a function of altitude but observed 
from the NASA shuttle. 
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Figure 10, Model used for predicted uncertainty calculations of DAS ozone concentration measure 
ments from the shuttle. 
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Figure 11. Predicted uncertainties as calculated for ozone concentration measurements per- 
formed from the NASA shuttle. 
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Figure 12. Ozone absorption coefficients vs wavelength 
reported by Griggs. 
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Figure 13. Predicted uncertainties as calculated for ozone column content measurements to 
various altitudes from the NASA shuttle. 
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